Pseudogap temperature as a Widom line in doped Mott insulators 
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The pseudogap refers to an enigmatic state of matter with unusual physical properties found below 
a characteristic temperature T* in hole-doped high-temperature superconductors. Determining T* 
is critical for understanding this state. Here we study the simplest model of correlated electron 
systems, the Hubbard model, with cluster dynamical mean-field theory to find out whether the 
pseudogap can occur solely because of strong coupling physics and short nonlocal correlations. 
We find that the pseudogap characteristic temperature T* is a sharp crossover between different 
dynamical regimes along a line of thermodynamic anomalies that appears above a first-order phase 
transition, the Widom line. The Widom line emanating from the critical endpoint of a first-order 
transition is thus the organizing principle for the pseudogap phase diagram of the cuprates. No 
additional broken symmetry is necessary to explain the phenomenon. Broken symmetry states 
appear in the pseudogap and not the other way around. 



Introduction 

The phase diagram of hole-doped high-tcmpcrature su- 
perconductors remains puzzling. A state of matter with 
unusual physical properties, dubbed "the pseudogap", is 
found below a characteristic temperature T* in a doped 
Mott insulator. Since the superconducting state is born 
out of the pseudogap over much of the phase diagram, 
the nature of the pseudogap is a fundamental issue in 
the field and it is under intense theoretical [1] and exper- 
imental [2-4] scrutiny. 

A pseudogap can occur because of disorder-broadened 
long-range ordered phases of Ising type [o] or because 
of fluctuating precursors to a long-range ordered phase 
that would appear only at T = because of the 
Mermin- Wagner theorem [(>] . The phenomenology of the 
pseudogap in the less strongly-coupled electron-doped 
cuprates [7, N] differs from that of hole-doped cuprates 
and is consistent with fluctuating precursors [!)]. The 
above two generic mechanisms have in common that they 
can occur even when interactions between electrons are 
not strong enough to lead to a Mott insulator at half- 
filling. Various broken symmetry states, such as spin- 
charge density wave, have been linked to the pseudo- 
gap in some cuprate families [2-4, 9] but not in all of 
them. Yet the pseudogap is a generic feature of all hole- 
doped cuprates. It is also possible that broken symme- 
tries are different phenomena, and in particular only a 
consequence, and not the origin, of the pseudogap, as 
suggested recently [] ()]. 

On the other hand, the Mott phenomenon, a blocking 
of charge transport because of strong electronic repul- 
sion, is ubiquitous in hole-doped cuprates [11] and it is 
appropriate to ask whether there is a third generic mech- 
anism for the pseudogap that is associated purely with 
Mott physics in two dimensions [7]. This mechanism ap- 



pears only if the Coulomb repulsion is strong enough [7] 
to turn the system into a Mott insulator at zero doping, 
as observed in hole-doped cuprates. No broken transla- 
tional or rotational symmetry is needed, although broken 
symmetry may also occur in certain cases [5]. In previ- 
ous work [12, 13], we identified a first-order transition 
at finite doping between two different metals. Here we 
show that the characteristic temperature T* is an unex- 
pected example of a phenomenon observed in fluids [14], 
namely a sharp crossover between different dynamical 
regimes along a line of thermodynamic anomalies that ap- 
pears above that first-order phase transition, the Widom 
line [14]. 

Results 

Widom line in doped Mott insulator. To investigate the 
formation of the pseudogap upon doping a Mott insu- 
lator, we study the competition between nearest neigh- 
bors hopping t and screened Coulomb repulsion U em- 
bodied in the two-dimensional Hubbard model. For the 
cuprates, the clear momentum dependence of the self- 
energy observed in photoemission (ARPES) forces one 
to use cluster extensions [15, IG] of dynamical mean- 
field theory [17]. It is known that four sites [LS] suf- 
fices to reproduce qualitatively the ARPES spectrum. 
Larger cluster sizes [19] will improve momentum reso- 
lutions and change quantitative details but should not 
remove first-order transitions, as has been verified at 
half- filling [20]. As in any cluster mean- field theory (dy- 
namical or not), we can study the normal-state phase 
of the model by suppressing long-range magnetic, super- 
conducting or other types of order while retaining the 
short-range correlations. Previous works based on this 
theoretical framework have shown that a pseudogap ap- 
pears close to a Mott insulator, as a result of short-range 
correlations [7, 18, 21, 22]. But it is only with recent the- 
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orctical and computational advances [23] that detailed 
results in the low temperature region of the phase dia- 
gram can be obtained. 

Recently, we used cellular dynamical mean-field theory 
on a 2 X 2 plaquette to demonstrate that a first-order tran- 
sition inhabits the finite doping part of the normal-state 
phase diagram of the model [f2, 13]. That first-order 
transition occurs along a coexistence line between two 
types of metals and ends at a critical point (Tp,/Xp) (see 
T — /i phase diagram in Fig. la). Here we show that the 
pseudogap is the low-doping phase whose characteristic 
temperature T* can be interpreted as the Widom line. 

The Widom line is defined as the line where the max- 
ima of different thermodynamic response functions touch 
each other asymptotically as one approaches the critical 
point [14]. The maxima become more pronounced on ap- 
proaching Tp, diverging at Tp. Fig. lb shows the chem- 
ical potential dependence of the charge compressibility 
K = l/n?{dn/d^)T for several temperatures above Tp. 
Far above Tp, k develops a maximum that increases and 
moves to higher doping with decreasing temperatures, in- 
dicating that the charge compressibility diverges at the 
critical point. The loci of the charge compressibility max- 
ima, maxj^K, are shown in Fig. la and give an estimate 
of the Widom line. 

Crossing of the Widom line involves drastic changes in 
the dynamics of the system, as indicated in previous in- 
vestigations on the phase diagram of fluids [14, 24]. Sim- 
ilarly, here we show that the pseudogap, as seen in the 
single-particle density of states and in the zero frequency 
spin susceptibility, arises at high temperature from cross- 
ing the Widom line that radiates out of the critical point. 

Identification of T* on the basis of the local density 
of states. First we study the development of the pseu- 
dogap in the local density of states A{uj), which can 
be accessed by tunnelling or photoemission spectroscopy, 
along paths at constant temperature or at constant dop- 
ing. Figs. 2(a,c) show the evolution of the density of 
states with doping at a fixed temperature above and 
below Tp respectively. In the Mott insulating state, 
at (5 = 0, the density of states consists of lower and 
upper Hubbard bands separated by a correlation gap. 
Upon hole doping, there is a dramatic transfer of spec- 
tral weight from high to low frequency, as a consequence 
of strong electronic correlations. The low frequency 
part of A{u}) develops a pseudogap, i.e. a depression 
in spectral weight, between a peak just below the Fermi 
level and a peak above the Fermi level. The pseudo- 
gap exhibits a two-peak profile that is highly asymmet- 
ric [22], reflecting the large particle-hole asymmetry ob- 
served experimentally [25]. Upon increasing doping, the 
particle-hole asymmetry decreases, the spectral weight 
inside the pseudogap gradually fills in, and the distance 
between the two peaks slightly decreases. Finally, at 
a temperature-dependent doping, the pseudogap disap- 
pears and a rather broad peak appears in the density of 




FIG. 1: Phase diagram and Widom line, (a), temper- 
ature T - chemical potential /i phase diagram of the two- 
dimensional Hubbard model in the normal-state obtained by 
cellular dynamical mean-field theory. We take units where 
a = t = R = fcs = 1. Here, f/ = 6.2 > [/mit ~ 5.95 
necessary to create a Mott insulator at half filling (/i — 0). 
The phase transition between two metals, the pseudogap 
and the correlated Fermi liquid, occurs along the coexis- 
tence line (dashed green line) that ends at the critical point 
{fip,Tp) = (—0.707,1/65) (orange circle). The coexistence 
line is bounded by the spinodals (blue lines with squares). 
The phase boundaries are obtained from scans at constant T 
monitoring the behavior of the occupation n versus jj, [12, 1^5]. 
The onset of the Mott insulator (black line with squares) is 
defined by a plateau in the occupation at n = 1. The com- 
pressible, metallic, regions have n < 1. Below Tp, the particle 
occupation! n changes discontinuously at the spinodal lines, 
indicating the first-order transition. The green coexistence 
line and the extrapolations to T = are a guide to the eye. 
The Widom line, i.e. the line where the maxima of differ- 
ent response functions converge, extends the coexistence line 
above Tp. We estimate this line by the loci of charge com- 
pressibility maxima max^K, where k = {dn/d^)T (red 
line with triangles), (b) Semi- logarithmic plot of the charge 
compressibility k, versus chemical potential /i for several tem- 
peratures above Tp, obtained by numerical derivative of the 
filling with respect to the chemical potential. The loci of max- 
ima of K are plotted in the T — fj, plane in (a). The value at 
the maximum increases as T decreases, indicating the diver- 
gence of K at the critical point. This is analog to the diverging 
charge compressibility at the Mott critical point [.i-i] or at the 
liquid-gas trans! ition. 



states which narrows with increased doping. The change 
from pseudogap to correlated Fermi liquid behavior oc- 
curs cither by a first-order transition when T < Tp, with 
a discontinuous change in A(lu), or by a crossover when 
T > Tp. We obtain T*{S) from the inflection point in 
A{uj = 0){6) at finite doping (see Figs. 2(b,d)). 

Figs. 2(e,g) show the temperature evolution of the den- 
sity of states at constant doping. In Fig. 2(g) above 6p 
(below fj.p), there is no pseudogap, only a peak around 
cj = that broadens with increasing temperature [13, 22]. 
By contrast, in Fig. 2(e) below 6p (above fip), the effect of 
increasing the temperature is to gradually fill the pseudo- 
gap, without decreasing the peak-to-peak distance. We 
identify the disappearance of the pseudogap in the spec- 
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FIG. 2: r* from density of states. (a,c,e,g), A{uj) at U — 
6.2 obtained from maximum-entropy analytic continuation of 
continuous-time quantum Monte Carlo data. (a,c), A{lu) for 
different dopings at fixed temperature, (a) T = 1/25 > Tp 
and (c ) T = 1/100 < Tp respectively. In the latter case, as- 
sociated with the first-order transition at finite doping there 
are hysteresis effects in the T — ^ plane, where the transi- 
tion appears as a coexistence interval between two metallic 
phases. Orange and magenta lines are A{u}) obtained in the 
coexistence region: they have the same chemical potential 
fj, = —0.725, but different filling and different low-energy spec- 
trum, one with a pseudogap and the other one with a rather 
broad peak at a; = 0. (e,g), density of states at fixed doping 
for various temperatures, (e) is for 5 = 0.02 < 5p, where a 
pseudogap opens up with decreasing temperature, (g) is for 
S — 0.08 > 5p where a pseudogap is absent and a narrow 
peak develops as T decreases. (b,d,f,h), density of states at 
the Fermi level A{(jj = 0) as a function of doping or temper- 
ature. Here, data are obtained from the extrapolated value 
of the imaginary part of the local cluster Green's function 
— l/7rImG(a;n 0), which does not require analytical contin- 
uation. The inflection point of these curves as a function of fi 
or T, indicated by a red circle, is our estimate of T*. On the 
upper horizontal axis we convert into physical units by using 
t = 0.35eV. 



trum by the inflection in A{u! — 0)(r) (see Fig. 2(f) and 
Supplementary Fig. SI). We will discuss later the evolu- 
tion with doping of T*, but the absence of pseudogap for 
(5 > already allows us to conclude that T* appears as 
a bridge between the Mott insulator and the first-order 
transition at finite doping. 

Pseudogap versus Mott insulator. Next we discuss the 
relationship between the pseudogap and the Mott in- 
sulator. The pseudogap is linked to Mott physics be- 



cause it opens up only above the critical U for the Mott 
transition, in the metal near the Mott insulator. A 
pseudogap can occur below this threshold due to long- 
wavelength antiferromagnetic fluctuations, but that is 
different physics [G] occurring at a different energy scale. 
Here, only short-range spin correlations are involved, as 
observed experimentally [20] in YBa2Cu408 at the pseu- 
dogap temperature. The first-order transition at finite 
doping, which is the terminus of the pseudogap phase, 
is linked as well to Mott physics because in the {U, T, /i) 
phase diagram, it emerges out of the Mott endpoint at 
half-filling, and progressively moves away from half-filling 
with increasing U [12, 13]. 

The pseudogap inherits many features from the parent 
Mott insulator: in both phases the electrons are bound 
into short-range singlets because of the superexchange 
mechanism, reminiscent of the resonating valence bond 
state [27]. Fig. 3a shows the temperature evolution of 
the probability obtained from the largest diagonal ele- 
ments of the reduced density matrix on a 2 x 2 clus- 
ter [2(S], which provide direct access to spin correlations 
(see also Supplementary Fig. S2). Below 6p, upon de- 
creasing temperature, the singlet is gaining weight at the 
expense of the triplet indicating a reduction of spin fluc- 
tuations. Because of singlet formation, the spin suscepti- 
bility xiT) drops below a characteristic temperature, as 
shown in Fig. 3b and as found in experiments [29]. The 
inflection point in xC^) defines a T* that moves to lower 
temperatures as the doping increases and approaches Tp 
as 5 Sp. 

Despite the magnetic behavior similar to the Mott in- 
sulator phase, the pseudogap phase is a new state of mat- 
ter. At T = it appears to be separated from the Mott 
insulator by a second-order transition [13]. With decreas- 
ing temperature, neither the value of the density of states 
at the Fermi level (Fig. 2(f)), nor the spin susceptibility 
(Fig. 3b), extrapolate to zero. The T* extrapolated to 
(5 = is not related to the opening of the Mott gap 
(Fig. 4). Finally, the peak to peak distance for the pseu- 
dogap does not extrapolate, as 5 — ?■ 0, to the Mott gap, 
just as is observed in experiment [30]. 

T* as the Widom line. Now we move to the relation- 
ship between the pseudogap and the Widom line ema- 
nating from the critical point at finite doping. Fig. 4 
shows in the T — 6 plane the doping evolution of the var- 
ious T*, identified above as infiection points in A{uj = 0) 
along constant T or constant 6 paths, and as inflection 
points in x(^) s-^d Probsingiot(T') at constant S. The 
different T* lines for the pseudogap move closer to each 
other along the Widom line as we approach the criti- 
cal endpoint. The interrelation between T* and Widom 
line is our main flnding. Therefore our work shows that 
the dynamic crossover associated with the buildup of the 
pseudogap is concomitant with a crossover in the ther- 
modynamic quantities, as observed in supercritical flu- 
ids. The organizing principle of these phenomena is the 
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FIG. 3: T* from spin susceptibility and plaquette 
eigenstates. (a), probability of the following plaquette 
eigenstates as a function of temperature, for several values 
of doping: the singlet \N = 4,5 = Q,K = (0,0)) and 
the triplet \N = A, S = 1, K = [11,11)) (squares and trian- 
gles respectively), where A*', S, K are the number of elec- 
trons, the total spin and the cluster momentum of the pla- 
quette eigenstate. (b), the zero-frequency spin susceptibility 
xiT) = {Sz{r)Sz{0))dT as a function of temperature for 
several values of doping. Sz is the projection of the total spin 
of the plaquette along the z direction. The inflection point of 
these curves as a function of T, marked by a solid symbol, is 
our estimate of T* . Data are for U = 6.2. 



Widom line. One can thus interpret T* as the Widom 
line, or, equivalent ly, consider the Widom line as a ther- 
modynamic signal of T*. Our results suggests that all 
indicators (both thermodynamic and dynamic) of the 
pseudogap temperature scale T* should approach each 
other with increasing doping, joining a critical endpoint 
of a first-order transition, which thereby appears as the 
source of anomalous behavior. In this view, pseudogap 
and strongly correlated Fermi liquid are separated from 
each other at low temperature by a first-order transition 
and are thus two distinct states of matter, just as liquid 
and gas are two distinct states, or phases. 
Discussion 

Common theories to explain the pseudogap phase in- 
clude the presence of rotational and/or spatial broken- 
symmetry phases as an essential ingredient. By contrast, 
in our approach the pseudogap is a consequence of large 
screened Coulomb repulsion that leads to strong singlet 
correlations in two dimensions reminiscent of resonating 
valence bond physics [27]. Competing phases are not nec- 
essary to obtain a pseudogap. The pseudogap phase can 
however be unstable to such phases. We therefore pro- 
vide a new and generic mechanism for the pseudogap in 
doped Mott insulators, according to which the pseudo- 
gap state is a new state of matter, whose characteristic 
temperature T* corresponds to the Widom line arising 
above a first-order transition. 

The Mott transition is often masked by broken- 
symmetry states. Similarly, our finite-doping transition 
is masked by the superconducting phase for instance [31]. 
Nevertheless the rapid crossover between pseudogap and 
metallic phases observed above the broken-symmetry 
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FIG. 4: Pseudogap temperature T* along the Widom 

line. T — 5 phase diagram for U = 6.2 (see Supplementary 
Fig. S3 for the T — jj, phase diagram). The pseudogap tem- 
perature scale T* is computed from the inflection point along 
paths at constant 5 of A(a; — 0)(line with green circles), x(.T) 
(magenta circles) and Probsingiet (r) (cyan circles). The T* 
line appears as soon as the Mott insulator is doped and joins 
the critical endpoint {Tp,5p) (orange circle) of a first-order 
transition away from half-filling, closely following the Widom 
line (red line with triangles) that emanates from the critical 
endpoint. The critical point moves to large doping and low 
temperatures with increasing U [12, 13]. The hatched region 
corresponds to the instability region bounded by the spin- 
odals. Extrapolations to T = are a guide to the eye. On 
the right vertical axis we convert into physical units by using 
t = 0.35eF. 



states is accessible and can be controlled by the Widom 
line. Such rapid change in dynamics is a hallmark of 
the Widom line [14, 24] and it is consistent with the 
strong-coupling nature and the observed phenomenology 
of the pseudogap in the vicinity of T* in hole-doped high- 
temperature superconductors. 

From a broader perspective, our work brings the con- 
ceptual framework of the Widom line, recently developed 
in the context of fluids [14, 24], to a completely different 
state of matter, the electronic fluid, suggesting its unex- 
pected generality. We recall the strong impact that re- 
sulted from bringing in the field of electronic properties 
of solids the well known concepts of smectic and nematic 
order developed earlier in the field of liquid crystals. It is 
tempting to argue that the same fate awaits the Widom 
line. 

Methods 

Our results are based on the cellular dynamical mean- 
field theory (CDMFT) [15, IG] solution of the two- 
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dimensional Hubbard model on the square lattice, 

ij(T i ^ ^ ^ ^ i 

(1) 

where c^o- and operators destroy and create electrons 
on site i with spin a, and = cf^Ci„. t is the hopping 
amplitude between nearest neighbors, U the energy cost 
of double occupation at each site and ^ the chemical po- 
tential. CDMFT isolates a cluster of lattice sites, here 
a 2 X 2 plaquette, and replaces the missing lattice en- 
vironment by a bath of non-interacting electrons which 
is self-consistently determined. The cluster in a bath 
problem is solved by a continuous-time Quantum Monte 
Carlo summation of all diagrams obtained by expand- 
ing the partition function in powers of the hybridization 
between bath and cluster [23, 28]. The size of the plaque- 
tte is large enough to be consistent with the experimental 
observation [2G] that at T*, in hole-doped cuprates, the 
antiferromagnetic correlation length is one or two lattice 
spacings. 

The value of Coulomb interaction U = &.2t is larger 
than the critical threshold J7mit ~ 5.95t necessary to ob- 
tain a Mott insulator at half-filling [12, 13] and is chosen 
such that the pseudogap critical temperature Tp{U) is 
accessible with our method. We carry out simulations at 
constant temperature for several values of fi and at con- 
stant doping 6 ~ 1 ^ n for several temperatures. Critical 
slowing down is a widespread and standard signal that 
the system is approaching a critical threshold [32] and 
appears in our simulations along the Widom line close to 
the critical point. To obtain reliable results, we perform 
up to 10^ Monte Carlo sweeps, averaged over 64 proces- 
sors, and hundreds of CDMFT iterations. The typical 
error on n is of order 10^^. 



[1] Norman, M. R., Pines, D. & Kallin, C. The pseudogap: 
friend or foe of high Tc? Adv. Phys. 54, 715 (2005). 

[2] Daou, R. et al. Broken rotational symmetry in the pseu- 
dogap phase of a high-Tc superconductor. Nature 463, 
519-522 (2010). 

[3] Lawler, M. J. et al. Intra-unit-cell electronic nematicity 
of the high-Tc copper-oxide pseudogap states. Nature 
466, 347-351 (2010). 

[4] He, R.-H. et al. From a Single-Band Metal to a High- 
Temperature Superconductor via Two Thermal Phase 
Transitions. Science 331, 1579-1583 (2011). 

[5] Fradkin, E., Kivelson, S. A., Lawler, M. J., Eisenstein, 
J. P. & Mackenzie, A. P. Nematic Fermi Fluids in Con- 
densed Matter Physics. Annual Review of Condensed 
Matter Physics 1, 153 -178 (2010). 

[6] Vilk, Y. M. & Tremblay, A.-M. S. Non-perturbative 
many-body approach to the hubbard model and single- 
particle pseudogap. J. Phys I (France) 7, 1309 - 1368 
(1997). 



[7] Senechal, D. & Tremblay, A.-M. S. Hot spots and pseudo- 
gaps for hole- and electron-doped high-temperature su- 
perconductors. Phys. Rev. Lett. 92, 126401 (2004). 
[8] Weber, C, Haule, K. & KotUar, G. Strength of cor- 
relations in electron- and hole-doped cuprates. Nature 
Physics 6, 574-578 (2010). 
[9] Motoyama, E. M. et al. Spin correlations in the 
electron-doped high-transition-temperature supercon- 
ductor NCCO. Nature 445, 186 (2007). 

[10] Parker, C. V. et al. Fluctuating stripes at the on- 
set of the pseudogap in the high-Tc superconductor 
BizSraCaCuaOg+x. Nature 468, 677-680 (2010). 

[11] Lee, P. A., Nagaosa, N. & Wen, X.-G. Doping a Mott in- 
sulator: Physics of high-temperature superconductivity. 
Rev. Mod. Phys. 78, 17 (2006). 

[12] Sordi, G., Haule, K. & Tremblay, A.-M. S. Finite 
Doping Signatures of the Mott Transition in the Two- 
Dimensional Hubbard Model. Phys. Rev. Lett. 104, 
226402 (2010). 

[13] Sordi, G., Haule, K. & Tremblay, A.-M. S. Mott physics 
and first-order transition between two metals in the 
normal-state phase diagram of the two-dimensional Hub- 
bard model. Phys. Rev. B 84, 075161 (2011). 

[14] Xu, L. et al. Relation between the Widom line and the 
dynamic crossover in systems with a liquid liquid phase 
transition. Proceedings of the National Academy of Sci- 
ences of the United States of America 102, 16558-16562 

(2005) . 

[15] Kotliar, G. et al. Electronic structure calculations with 
dynamical mean-field theory. Rev. Mod. Phys. 78, 865 

(2006) . 

[16] Maier, T., Jarrell, M., Pruschke, T. & Hettler, M. H. 
Quantum cluster theories. Rev. Mod. Phys. 77, 1027- 
1080 (2005). 

[17] Georges, A., Kotliar, G., Krauth, W. & Rozenberg, 
M. J. Dynamical mean-field theory of strongly correlated 
fermion systems and the limit of infinite dimensions. Rev. 
Mod. Phys. 68, 13 (1996). 

[18] Kyung, B. et al. Pseudogap induced by short-range spin 
correlations in a doped Mott insulator. Phys. Rev. B 73, 
165114 (2006). 

[19] Gull, E., Ferrero, M., ParcoUet, O., Georges, A. & Mib 
lis, A. J. Momentum-space anisotropy and pseudogaps: 
A comparative cluster dynamical mean-field analysis of 
the doping-driven metal-insulator transition in the two- 
dimensional Hubbard model. Phys. Rev. B 82, 155101 

(2010) . 

[20] Park, H., Haule, K. & Kotliar, G. Cluster dynamical 
mean field theory of the Mott transition. Phys. Rev. 
Lett. 101, 186403 (2008). 

[21] Macridin, A., JarreU, M., Maier, T., Kent, P. R. C. & 
D'Azevedo, E. Pseudogap and Antiferromagnetic Cor- 
relations in the Hubbard Model. Phys. Rev. Lett. 97, 
036401 (2006). 

[22] Haule, K. & Kotliar, G. Strongly correlated superconduc- 
tivity: A plaquette dynamical mean-field theory study. 
Phys. Rev. B 76, 104509 (2007). 

[23] Gull, E. et al. Continuous-time monte carlo methods for 
quantum impurity models. Rev. Mod. Phys. 83, 349-404 

(2011) . 

[24] Simeoni, G. G. et al. The Widom line as the crossover 
between liquid-like and gas-like behaviour in supercritical 
fiuids. Nature Physics 6, 503-507 (2010). 

[25] Kohsaka, Y. et al. An intrinsic bond-centered elec- 



tronic glass with unidirectional domains in underdoped 
cuprates. Science 315, 1380-1385 (2007). 

[26] Curro, N. J., Imai, T., Slichter, C. P. & Dabrowski, B. 
High-temperature ^''Cu(2) nuclear quadrupole and mag- 
netic resonance measurements of YBa2Cu408 . Phys. 
Rev. B 56, 877-885 (1997). 

[27] Anderson, P. W. The resonating valence bond state in 
La2Cu04 and superconductivity. Science 235, 1196 - 
1198 (1987). 

[28] Haule, K. Quantum Monte Carlo impurity solver for clus- 
ter dynamical mean-field theory and electronic structure 
calculations with adjustable cluster base. Phys. Rev. B 
75, 155113 (2007). 

[29] AUoul, H., Ohno, T. & Mendels, P. ^^Y NMR evidence 
for a fermi-liquid behavior in Y Ba2Cus,0(i+x. Phys. Rev. 
Lett. 63, 1700-1703 (1989). 

[30] Hiifner, S., Hossain, M. A., Damascelli, A. & Sawatzky, 
G. A. Two gaps make a high-temperature superconduc- 
tor? Reports on Progress m Physics 71, 062501 (2008). 

[31] Sordi, G., Semon, P., Haule, K. & Tremblay, A.- 
M. S. Strong Coupling Superconductivity, Pseudogap, 
and Mott Transition. Phys. Rev. Lett. 108, 216401 
(2012). 

[32] Sheffer, M. et al. Early-warning signals for critical tran- 
sitions. Nature 461, 53-59 (2009). 
[33] Kotliar, G., Murthy, S. & Rozenberg, M. J. Compressibil- 
ity Divergence and the Finite Temperature Mott Transi- 
tion. Phys. Rev. Lett. 89, 046401 (2002). 
Acknowledgments We thank E. Kats and L. Taille- 
fer for discussions. This work was partially supported 
by FQRNT, by the Tier I Canada Research Chair Pro- 
gram (A.-M.S.T.), and by NSF DMR-0746395 (K.H.). 
Computational resources were provided by CFI, MELS, 
the RQCHP, Calcul Quebec and Compute Canada. A.- 
M.S.T is grateful to the Harvard Physics Department 
for support and P.S. for hospitality during the writing 
of this work. Partial support was also provided by the 
MIT-Harvard Center for Ultracold Atoms. 

Author contributions G.S. conceived the project 
and carried the data analysis. P.S. and K.H. wrote the 
main codes. G.S. and A.-M.S.T. wrote the paper and 
all authors discussed it. A.-M.S.T. supervised the entire 
project. 

Competing financial interests The authors declare 
no competing financial interests. 



7 



Pseudogap temperature as a Widom line in doped 
Mott insulators: supplementary information 

G. Sordi, P. Semon, K. Haule, A.M.-S. Tremblay 



QQ max dA({D=0)/d^ 
Q£] max^dA((D^/dT 




FIG. SI: (a,b), density of states at the Fermi level A{io = 0) 
as a function of (a) chemical potential fi or (b) temperature 
T. Data are obtained from the extrapolated value of the local 
cluster Green's function l/nlmG{!jj„ — >■ 0). (c ), pseudogap 
temperature scale T* in the T — S phase diagram. We evaluate 
T* by the inflection point in A{ijj = 0)(/x) (line with black 
circles) and A{uj = 0)(T) (line with red squares). These two 
estimates of T* are consistent and terminate at the critical 
point (orange circle) of the first-order transition. At larger 
value of U [12, 13], one would obtain a larger value of 5p, 
but for large 5p the corresponding Tp becomes inaccessible 
numerically. 
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FIG. S3: T—jj, phase diagram. With respect to Fig. 4, only the 
X-axis has changed. All indicators of the pseudogap tempera- 
ture scale T* follow the Widom line (line with red triangles) 
and terminate at the critical point (orange circle). 
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FIG. S2: Probability of the following 2x2 plaquette eigen- 
states as a function of temperature, for several values of dop- 
ing: the singlet \N = i, S ^ 0, K ^ (0,0)) and the doublet 
|A'^ = 3, S = 1/2, K = (7r,0)) (squares and circles respec- 
tively), where A'^, S, K are the number of electrons, the total 
spin and the cluster momentum of the plaquette eigenstate. 
Although the spin doublet is more probable than the triplet 
over much of the temperature range, its probability is essen- 
tially temperature independent. 



